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ABSTRACT 

Linear polarization at the level of ^ 1 - 3% has by now been measured in several GRB afterglows. Whereas 
the degree of polarization, P, was found to vary in some sources, the position angle, 9p, was roughly constant 
in all cases. Until now, the polarization has been commonly attributed to synchrotron radiation from a jet with 
a tangled magnetic field that is viewed somewhat off axis. However, this model predicts either a peak in P or a 
90° change in 9p around the "jet break" time in the lightcurve, for which there has so far been no observational 
confirmation. We propose an alternative interpretation, wherein the polarization is attributed, at least in part, 
to a large-scale, ordered magnetic field in the ambient medium. The ordered component may dominate the 
polarization even if the total emissivity is dominated by a tangled field generated by postshock turbulence. In 
this picture, Oj, is roughly constant because of the uniformity of the field, whereas P varies as a result of changes 
in the ratio of the ordered-to-random mean-squared field amplitudes. We point out that variable afterglow light 
curves should be accompanied by a variable polarization. The radiation from the original ejecta, which includes 
the prompt 7-ray emission and the emission from the reverse shock (the 'optical flash' and 'radio flare'), could 
potentially exhibit a high degree of polarization (up to ^ 60%) induced by an ordered transverse magnetic field 
advected from the central source. 

Subject headings: gamma rays; bursts — polarization — radiation mechanisms; nonthermal — shock waves 



1. INTRODUCTION 

The first detection of polarization in an optical afterglow 
was in GRB 9905 10, where a degree of polarization P= 1 .7 ± 
0.2% (1.6 ±0.2%) was measured at fobs « 18 hr (21 hr) after 
the burst. ^ Since then, P ^ 1-3% has been detected in a 
few additional afterglows,^ some of which showed a temporal 
variation in P (Rol et al. 2000; Barth et al. 2003), but typically 
the position angle (PA) 6p showed little or no change. A few 
other afterglows produced only upper limits, P < 2-5%. 

The polarization is attributed to synchrotron emission be- 
hind a shock wave. It thus depends on the local magnetic 
field configuration, which determines the polarization at each 
point of the afterglow image, and on the global geometry of 
the shock, which determines how the polarization is averaged 
over the (unresolved) image. We make a distinction between 
magnetic field configurations that are axially symmetric about 
the normal «sh to the shock surface and those that are not.^ 
The first category gives no net polarization for a spherical 
flow and is assumed in most previous works (Sari 1999; Ghis- 
elhni & Lazzati 1999; Medvedev & Loeb 1999; Granot et al. 
2002; Rossi et al. 2002). These models take the field to be 
completely random in the plane of the shock, and the polar- 
ization is usually attributed to a jet viewed somewhat off axis. 
For a structured jet, P has one peak near the jet break time 
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There is one exception, P = 9.9 ± 1.3%, measured in GRB 020405 at 
fobs = 1.3 days (Bersier et al. 2003). Significantly lower values (P Ri 1.5 — 
2%) were measured in this afterglow at t^bs = 1-2,2.2, and 3.3 days by other 
groups, with a similar 9p . If real, this behavior has no simple explanation in 
any of the existing models. 

^ In this picture, the field is tangled over very small scales and possesses 
this symmetry when averaged over regions of angular size <C 1 /7, where 7 
is the Lorentz factor of the shocked fluid. 



tj (when 1 /7 increases to 6q, the initial jet opening half- 
angle), whereas for a uniform jet P has 2 (or 3) peaks near tj, 
with P passing through zero and 6p changing by 90° between 
the peaks. 

The second category can produce net polarization even for 
a spherical flow. One example is the "patchy coherent field" 
model of Gruzinov & Waxman (1999), where the observed 
region consists of 50 mutually incoherent patches of an- 
gular size < 1 /7, within each of which the field is fully or- 
dered.'' This model predicts P - P^^^/N^'^ 10% (where 
^max ^ 60-70% is the maximum P of local synchrotron emis- 
sion in a uniform magnetic field) and simultaneous (random) 
variability in P and dp on timescales Afobs fobs ■ The idea 
behind this model is that ^ 1 /7 is the angular size of causally 
connected regions, and for a magnetic field that is generated 
in the shock itself this is the largest scale over which the field 
can be coherent. 

However, if the magnetic field were ordered on an angu- 
lar scale I/7, then the resulting P could approach Pmax- 
Such a situation can be realized if an ordered field exists in 
the medium into which the shock propagates. For a typical 
interstellar medium, the postshock field would still be very 
weak (with the magnetic energy a fraction eg ^ 10"'" of the 
internal energy), but it would be higher {eg < 10""^; Biermann 
& Cassinelli 1993) if the shock expands into a magnetized 
wind of a progenitor star, and higher yet (eg ^ 0.01 -0.1) if it 
propagates into a pulsar-wind bubble (PWB), as expected in 
the supranova model (Konigl & Granot 2002). A strong or- 
dered field component is likely to exist in the original ejecta 
and could give rise to a high value of P (< P^ax) in both the 
prompt GRB and the reverse-shock emission. 

We calculate the polarization for a jet with a tangled mag- 
netic field in § 121 investigate the effects of adding an ordered 
field component in §|3] and discuss the results in §|4] 

* A similar model was used to study the linear polarization induced by 
microlensing of GRB afterglows (Loeb & Pema 1998). 
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2. POLARIZATION FROM A JET WITH A TANGLED MAGNETIC 

FIELD 

Synchrotron emission is generally partially linearly po- 
larized. In terms of the Stokes parameters: V = 0, Op = 
i arctan(t//0, and P = {Q^ + U^)^l^/I. As the Stokes param- 
eters are additive for incoherent emission, they can be cal- 
culated by summing over all the contributions from different 
fluid elements to the flux at a given observed time fobs- In 
practice, the flux is calculated by dividing fobs into bins of size 
(5fobs, and assigning to the appropriate time bin the contribu- 
tion to F„ from the emission of each 4-volume fluid element 

- ^_ (1+^)' l'(r,tmStobs-2\t-^-t,^,,\)d'x 
72(r,0[l-«-v(r,0/c] Jfobs 

(1) 

[q,i]-{T.^^^) E^^^|peos2i|' 

where f^, is the local emissivity, v is the fluid velocity, h is 
the direction to the observer, Q is the Heaviside step function 
and the summation in Eq. |2]is over r & f, for fixed fobs & «■ 

In this section, we consider only a random magnetic field 
Brnd that is tangled on angular scales <C 1 /7, with axial sym- 
metry w.r.t. Hsh- The field anisotropy is parameterized by 
b = 2{BJ^) / {B\), where B|| (B±) is the magnetic field com- 
ponent parallel (perpendicular) to fish- The local polarization 
of the emission from a given fluid element is^ 

Prnde')_ {{Bi)-{Bl)/2)sm'e' _ (t-Dsin'e' 
Pm-ax {Bj } sin^ 61' + ( 1 + cos2 6*') (Bi ) /2 2 + (Z? - 1 ) sin^ 6*' ' 

(3) 

(Gruzinov 1999; Sari 1999), where cos6'' = = h-h,^. 

For P > 0, the polarization is in the direction of n x Wsh- 

The magnetic field configuration behind the shock, and 
hence the value of b, cannot be easily deduced from first prin- 
ciples. It was suggested that small-scale postshock fields can 
be generated by a two-stream instability (e.g., Medvedev & 
Loeb 1999), which predicts b <^\. However, it is not clear 
whether the magnetic fields produced in this way survive in 
the bulk of the postshock flow (Gruzinov 1999, although see 
Frederiksen et al. 2003), or whether this is the dominant tan- 
gling mechanism. The relatively low observed values of P 
(< 3-4%) suggest that 0.5 < < 2 if the polarization is due 
to a jet with a shock-generated field. 

Turbulence in the postshock region (possibly induced by a 
microinstability) could amplify and isotropize the field, keep- 
ing b close to 1. As each fluid element moves downstream 
from the shock transition, it is sheared by the flow. For a 
Blandford-McKee (1976) self-similar blastwave solution with 
an ambient density pext « r"*^, the length of a fluid element in 
the directions parallel and perpendicular to Wsh scales with the 
self-similarity variable x (where x = 1 at the shock front and 
increases with distance behind the shock) as L|| oc ;^2(4-t) and 

L_L oc x*^. Therefore, the stretching of each fluid element 
in the radial direction is larger than in the tangential direc- 
tion. This would increase b while maintaining axial symme- 
try about «sh- The relevant value of b here is the average over 

^ For simplicity, we assume j'^, oc (B' sinx')'', where cosx' = ■ B' , with 
e = 2 (Sari 1999), although in reality e may be different. 
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Fig. 1. — Polarization lightcurves for a jetted GRB afterglow with a ran- 
dom magnetic field. The solid (dashed) lines are for h = (b = oo, with a 
minus sign); Pgo = f(0.6/Pmax). The jet parameters are the same as in Fig. 
3 of Granot et al. (2002). The lower right inset shows P* = sgn(l -b)P, 
normalized to its maximum value, for a viewing angle Oabs = 0.5 Bg and 
h = 0,0.5,2,00, making it easier to follow the effect of b on the shape of 
the lightcurve. The upper left inset shows a schematic diagram of the plane 
of the sky. The shaded region represents a jet with both a tangled and ordered 
field components. The projection of the ordered magnetic field on the plane 
of the sky, Bord, is at an angle S in the counterclockwise direction w.r.t. the 
direction from the line of sight (l.o.s.) to the jet axis. The polarization vector 
P is at an angle 6p, measured clockwise from the perpendicular to Bquj. 

the postshock region, weighted by the emissivity. If the turbu- 
lence only persists over a small distance (0 < x~ 1 ^ 1> but 
still the plasma skin depth) then, since most of the emission 
is from x a few, we may have I < b <a few.** 

Figure ^ shows the polarization lightcurves for b = and 
oo, based on the jet model of Kumar & Panaitescu (2000). 
For all viewing angles ^obs < Oq there are two peaks in P, a 
little before and after tj: P passes through zero in between 
these peaks as 9p changes by 90°. This result is similar to that 
of Ghisellini & Lazzati (1999), who did not consider lateral 
spreading of the jet, and differs from that of Sari (1999), who 
assumed 6'jet(fobs>0) =1/7' since the lateral spreading in the 
jet model that we use is smaller than the one used by Sari. 
The main distinction between the b < I and b > I cases is a 
90° difference in 9p, but this prediction can only be tested if 
one can independently determine the direction from the l.o.s. 
to the jet axis. This may in principle be done by measuring 
the direction of motion of the flux centroid (Sari 1999): for 
b < 1, P is perpendicular to (aligned with) n^. before (after) tj, 
whereas for b > I the situation is reversed. More generally, 
one would then be able to test if indeed P \\ He or P 1. he, as 
expected for a pure Bmd field, or if the angle between them 
is different, as would generally be the case if an ordered field 
component were also present. 

3. THE EFFECTS OF AN ORDERED MAGNETIC FIELD 

We now add an ordered magnetic field component Bo^i to 
the random field Bnid considered in §|2l In passing through the 
shock transition, the parallel component of the ambient mag- 
netic field Bext remains unchanged but the transverse compo- 
nent is amplified by a factor equal to the fluid compression 
ratio, which for 7 :§> 1 is 47. Thus typically B±^ ^ Z?|| behind 
the shock. For simplicity we assume that Bord lies in the plane 
of the shock and is fully ordered and that Bgxt is uniform, so 
that Bord is coherent over the entire shock. 

* On the other hand, if an intrinsically isotropic turbulence persists over a 
large portion of the emission region, this could reduce the effect of shearing 
on b, resulting in0<(6— 1)<C1. 
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Fig. 2. — Polaiization lightcurves for ordered + random magnetic field 
components, for b = and the same jet parameters as in Fig. As dp(5 = 
0,90°) = 0, we show only dp(t,5 = 45°). 
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Fig. 3. — Polarization lightcurves for a jet with 6q = 3°, 0obs - O.750o. 
h = Q.5, es tot = 0.05, ee =0.1, tj fafo = 14 hr, and Tj(t) given by eg o = 0-0015, 
a = 2, j(: = 6, g = 3, and5 = 3. 



It is most convenient to sum over the Stokes parameters as- 
sociated with Brnd and Boid separately and combine them at 
the end.^ The direction of polarization P^rd of the emission 
from the ordered component is perpendicular to its projection 
(Bold) on the plane of the sky; Pmd is either along the plane 
containing the l.o.s. and the jet symmetry axis (for P^d > 0) 
or perpendicular to that direction (for Pi„d < 0). The total po- 
larization and the PA are given by 
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l + r] 



= - arctan 
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where 77 = /oid/A-nd ~ (^ord) / (^rnd) ^^^^'^ °f obscrvcd 
intensities in the two components and 9p , 5 are measured as 
illustrated in the upper left inset of Figure^ 

For Z?ord we find that the PA as a function of the polar angle 
9 from the l.o.s. and the azimuthal angle (f) (measured from 
Boi-d) is given, in the relativistic (7^1) limit, by 9p = 4> + 
arctan(|^ cot0), where y = ("f9)^. We have Ii, = Il,{v/v')^, 
with /;;,, oc v'-'^{l-in' -B'^^ifyl^, where v/v' « and 



I+y 



are given by {U , Q)/IP„, 



i+sin 0. The Stokes parameters 
: / £/n/i.(sin26'p,cos26ip)/ / dVlU. 



For a spherical flow or a jet dX t <tj, when the edge of the 
jet is not visible, J dVl = J^^ dcf) dy, U = and 9p = 
Tr/2. For e = 2 and jmax = 1 we obtain -Q/IPmax = /(a) = 

- so Pord/i'max = /(") « 0.90-0.93 



(2^+°-I)(2+a)(3+a) 



27+Q_28+2a[23+a5-7+(23+°-l)a] 



for < a < 1.5. For a = ^, P„ax = -S^, e = 4^ (i.e. PLS 



2 . max y,+7/3' " 2 

G in Granot & Sari 2002), our analytic result corresponds to 
p = 3, for which Pord/^'max 0.92 and = « 0.69. 

For ji^ax » 1 , e = 2 we obtain /(a) = ^-j^^ ; for a = ^ and 
P = 3, Pord/Pmax = f » 0.86 and P = ^ « 0.64. The diference 
between the y^ax = 1 and y^ax ^ 1 results may be relevant to 
the prompt GRB (see § |^, where the tail of a pulse (which 
corresponds to y > 1) is predicted to be less polarized than its 
peak (y < 1). When the edge of the jet is visible, the limits 

' This is valid in the limit where the two components are associated with 
distinct fluid elements. Alternative schemes for combining Bjnd and Bord may 
produce a somewhat different polarization. 



of integration over dil change. As this causes relatively small 
modifications in 9p and Pord, we use the analytic expressions 
above for simplicity. 

Figure |2] depicts a sample of polarization lightcurves in 
which both 77 and b are taken to be independent of time. In 
this case the Boid-induced polarization is constant (in both P 
and 9p) throughout the afterglow. Interestingly, a similar po- 
larization signature could be produced by dust in our galaxy or 
in the GRB host galaxy. In the latter case, however, the polar- 
ization would likely be accompanied by absorption that would 
redden the spectrum: this could in principle make it possi- 
ble to estimate the level of the galactic dust contribution and 
thereby determine the fraction of such a constant-polarization 
component that is intrinsic to the source. Since Pord ~ Pmax 
whereas P,nd is typically much smaller, we find that for 77 = 1, 
and even for 77 = 1/3, the polarized intensity is still dominated 
by Bord (vPord > Prnd), with Bmd Only inducing relatively small 
fluctuations around the Z?oid-induced values of P and 9p. For 
6 = 45° the fluctuations in both P and 9p are very small in this 
parameter range. 

If Bold dominates the polarization, then, by equation (|4}, the 
time evolution of P follows that of the intensity ratio 77. The 
low measured values of P indicate that 77 ^ 1, so Bnid dom- 
inates the emissivity. To the extent that the random field is 
close to equipartition (e^.^d 1), 77 = eB.ord/cB.md ~ eB,ord- If 
the shock is radiative during its early evolution, then cooling- 
induced compression increases the emissivity-weighted e^ ord 
over its immediate postshock (adiabatic) value by a factor'" 
^ (1 -£(.)"'. The transition from fast to slow cooling, which 
occurs at f = to, could therefore reduce 77 and may contribute 
to the early decline of P observed in some sources. Dur- 
ing the subsequent, slow-cooling phase, e^^ord is essentially 
equal to the magnetization parameter of the ambient medium, 

= Bgx,/47rpextC^, so the evolution of P during that phase may 
reflect the radial behavior of this parameter: a is expected to 
be roughly constant for an ISM or a stellar wind but to in- 
crease with r inside a PWB (Konigl & Granot 2002). If the 
orientation of the ambient field also changed with radius then 
this would lead to a gradual variation in 9p. If one approxi- 

mates Bgxt oz r"' and pext « r , then e^.ord « f^'j^^ . We param- 
eterize the above effects by ejj^ord = eB,o(fobs/fo)^-F(fobs/fo), 

This assumes that the fraction of the internal energy just behind the 
shock transition that resides in relativistic electrons and pairs is radiated 
away (Granot & Konigl 2001). 
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where F(x) = 1 + arctan(^lnx) describes the ampHtude 
(g) and sharpness (0 of the change in eB,ord at fobs ^ h). We 
assume that e^.tot = eB,oi-d + eB,rnd = const, so ry = - I)"'- 

Figure |3l shows an example of the polarization lightcurves. 
The choice of parameters for 77(f) is motivated by GRB 
020813, in which 6p was again roughly constant with time 
whereas P first decreased (from P « 2% after ~ 6hr to P « 
0.6% after ^ 24 hr) and subsequently increased monotoni- 
cally, reaching P w 3.7% after 96hr. A sharp break in the 
lightcurve was observed after 14 hr (Covino et al. 2003b). 
It is seen that roughly equal contributions to the polarization 
from Z?ord and Bmd can provide a qualitative fit to the evolution 
of P and 9p in this source for S « 60° -90°. More generally, 
the polarization lightcurves can show a diverse behavior that 
varies as a function of 6 as well as of b and Oobs/So- So long 
as rjPoid > Prnd, the changes in 0p would be small whereas the 
variations in P could be significant, as found observationally. 

4. DISCUSSION 

The linear polarization in GRB afterglows may be largely 
due to an ordered magnetic field in the ambient medium, 
which gives rise to an ordered field component behind the af- 
terglow shock that is coherent over the entire emission region. 
This can result in a polarization position angle (PA), 9p, that 
is roughly constant in time as well as in a variable degree of 
polarization, P, as found in all afterglow observations to date 
(except one; see footnote [T3t . 

The magnetic field in the GRB ejecta is potentially much 
more ordered than in the shocked ambient medium behind the 
afterglow shock, reflecting the likely presence of a dynami- 
cally important, predominantly tansverse, large-scale field ad- 
vected from the source (e.g.. Spruit, Daigne, & Drenkhahn 
2001; Vlahakis & Konigl 2001). This could result in a large 
value of P [up to - (0.90-0.93)Pmax 60%] in the prompt 
7-ray emission" as well as in the 'optical flash' and 'radio 
flare', which are attributed to emission from the reverse shock. 
If the polarization from the reverse shock is indeed dominated 
by the ordered component, and if it is coherent over the whole 
ejecta, then 9p is not expected to vary significantly during 



the optical flash or between the optical flash and the radio 
flare. However, if the ordered magnetic field is coherent only 
in patches of angular size Og, then, so long as 7 > 1 /Og, we 
expect P Pmax, whereas after 7 drops below 1 /^^ we expect 
P ^ jOsPm-dx and variations in 6p on timescales Afobs ^ fobs 
on account of the averaging over ~ (i^bT^ mutually inco- 
herent patches within the observed region of angle 1 /7 about 
the line of sight. (This resembles the proposal by Gruzinov 
& Waxman 1999, except that here is envisioned to increase 
with time.) In the latter case P might be smaller and 9p would 
be different in the 'radio flare' (for which typically 7 < 10) 
than in the 'optical flash' (for which 7 > 100). 

Variability in the afterglow lightcurve, as reported in GRBs 
021004 and 030329, whether induced by a clumpy external 
medium or a patchy shell (Lazzati et al. 2002; Nakar, Piran, 
& Granot 2002), should give a different weight to emission 
from different parts of the afterglow image, thus breaking its 
symmetry and inducing polarization.'^ Therefore, we expect 
a highly variable lightcurve to be accompanied by variability 
in both P and 6*^.'^ 

Early polarization measurements, starting at fobs ^ tj, are 
crucial for distinguishing between our model and purely tan- 
gled jet field models, as the latter predict P(fobs ^ tj) ^ 
Pitohs ^ tj), whereas our model allows P(fobs ^ tj) ^ P(fobs ~ 
tj). In the latter models P is expected to peak, or else van- 
ish and reappear rotated by 90°, around tj. In contrast, in our 
model, if the polarization is dominated by an ordered mag- 
netic field, then the variations in the polarization around tj 
would be much less pronounced, with 9p exhibiting only a 
gradual variation and P never crossing zero. Our model pre- 
dicts a possible change in P around the transition time from 
fast to slow cooling, fo, where typically fo ~ 1 hr ('^ 1 day) for 
ISM-like (stellar wind-like) parameters (although it may vary 
considerably around these values). 

We thank P. Goldreich, R. Sari, A. Panaitescu, and E. Rossi 
for useful discussions. This research was supported in part by 
funds for natural sciences at the Institute for Advanced Study 
(JG) and by NASA ATP grant NAG5-12635 (AK). 



' ' After this paper was submitted, Coburn & Boggs (2003) reported a mea- 
surement of P = 80 ± 20% in the 7-ray emission of GRB 021206, which is 
naturally (and most likely; Granot 2003) produced in this way. 

If the density distribution is spherically symmetric, the symmetry would 
need to be broken by the outflow geometry — e.g., a jet observed off axis. 

After this paper was submitted, a change of 45° in 8p was reported in 



GRB 021004 between 9 and 16 hr (Rol et al. 2003). This cannot be explained 
by simple jet models (Sari 1999; Ghisellini & Lazzati 1999) but could natu- 
rally arise in conjunction with the variability in the lightcurve (which, in fact, 
peaked at about the same time). 
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